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We investigate the effective interaction of a pseudoscalar particle with two photons in an external 
electromagnetic field in the general case of all external particles being off the mass shell. This inter- 
action is used to study the radiative decay a — > 77 of the axion, a pseudoscalar particle associated 
with the Peccei-Quinn symmetry, and the crossing channel of photon splitting 7 ~~* ya. It is shown 
that the external field removes the suppression associated with the smallness of the axion mass and 
so leads to the strong catalysis of these processes. The field-induced axion emission by the photon 
is analyzed as one more possible source of energy losses by astrophysical objects. 
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I. INTRODUCTION 

For quite a long time there has been interest in various 
extensions of the standard model (SM) with a nonstan- 
dard set of Higgs particles. In the simplest version of the 
SM one needs only one Higgs doublet. However, there 
are reasons to believe that if scalar bosons exist at all, 
their number may be significant. The Higgs sector can 
be so rich that a situation becomes possible when one or 
several Goldstone degrees of freedom are not absorbed 
by the Higgs mechanism and appear as physical mass- 
less particles. Such Goldstone (or pseudo Goldstone) 
bosons, whose existence is related to the breakdown of 
exact (approximate) symmetries, is also of interest from 
\ the experimental point of view: they could be observed 
in low-energy experiments. Such particles include, e.g., 
axions lip], Majorons ||, arions H. The most widely 
discussed pseudo Goldstone boson, the axion, associated 
with the Peccei-Quinn (PQ) symmetry continues to be 
an attractive solution to the problem of CP conservation 
in strong interactions |). Although the original axion, 
associated with the spontaneous breakdown of the global 
PQ symmetry at the weak scale, f w , is excluded experi- 
mentally, modified versions of PQ models with their asso- 
ciated axions are of great interest. If the breaking scale 
of the PQ symmetry, f a , is much larger than the elec- 
troweak scale f a 3> f w , the resulting "invisible axion" is 
very weakly interacting (coupling constant ~ /a -1 ), very 
light (m a ~ fa 1 ) arL d very long lived. The axion lifetime 
in vacuum is gigantic: 



which leave a rather narrow window §-||Q: 

icr 5 cV<TO a < icr 2 ev, (2) 

where axions could exist and provide a significant frac- 
tion or all of the cosmic dark matter. A survey of vari- 
ous processes involving the production of weakly coupled 
particles and astrophysical methods for obtaining con- 
straints on the parameters of axion models is given in 
Ref. @. 

The important role which axions play in elementary 
particle physics stimulates constant interest in theoret- 
ical and experimental studies. Along with laboratory 
searches |l2,[l3|], stars are a unique laboratory to inves- 
tigate physical processes jyj. For example, the upper 
bound on m a is obtained from the requirement that stars 
not lose too much energy by axions. In studies of pro- 
cesses occurring inside astrophysical objects one has to 
take into account that dense matter substantially influ- 
ences particle properties. An intensive electromagnetic 
field also plays the role of an active medium modifying 
the dispersion relations of particles so that novel pro- 
cesses are not only opened kinematically, but become 
substantial as well. Photon splitting 7 — ► 77 Q] and 
the Cherenkov process v — > ^7 |jl5 16 are among them. 

Actually, both components of the active medium, 
plasma and a magnetic field, are presented in most as- 
trophysical objects. A situation is also possible when the 
magnetic component dominates. For example, in a su- 
pernova explosion or in a coalescence of neutron stars a 
region outside the neutrinosphcre of order of several tens 



r(a -> 2 7 ) - 6.3 x 10 48 s 



1(T 3 eV 



E„ 



lMeV 



(1) 



The allowed range for the axion mass is strongly con- 
strained by astrophysical and cosmological considerations 



* However in Ref. |10| a possibility to solve the CP problem 
of QCD within a grand unified theory model with a heavy 
axion m a < 1 TeV is considered. 
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of kilometers with a strong magnetic field and a rather 
rarefied plasma could exist. The possible existence of 
astrophysical objects with B ~ 10 15 — 10 17 G, signifi- 
cantly above the critical, Schwinger value B e = rnr je ~ 
4.41 x 10 13 G, was discussed both for toroidal and 
for poloidal || fields. 

In this paper we investigate a field-induced effective 
interaction of a pseudoscalar particle with two photons 
described by three-point loop diagrams. By now the de- 
tailed studies of the analogous three-point loop vertex 
of the photons' interaction in external electromagnetic 
fields are known only. The history of these investigations 
goes back to the pioneer paper by Adler and it is still 
in progress ]19|-|22;1. 

Note that amplitudes of processes in an external 
field depend not only on the kinematical invariants of 
type m? or p 2 , but also on field invariants |e 2 (FF)| 1 / 2 , 
|e 2 (FF)| 1 / 2 , |e 2 (pFFp)| 1 / 3 where m andp M are the mass 
and the four-momentum of a particle, is the external 
field tensor, and (pFFp) — p^F^ l/ F up p p . In the ultrarel- 
ativistic limit the field invariant \e 2 (pFFp)\ 1 / 3 can occur 
as the largest one. This is due to the fact that in the 
relativistic particle rest frame the field may turn out to 
be of order of the critical one or even higher, appearing 
very close to the constant crossed field (E _L B, E = £>), 
where (FF) = (FF) = 0. Thus, calculations in this field 
are the relativistic limit of calculations in an arbitrary 
relatively smooth field. 

In Sec. II we present the result of our calculations of 
the 077 vertex in the external crossed field in a general 
case when the external particles are off the mass shell. In 
Sec. Ill we discuss photon and axion dispersion relations 
in the crossed field and analyze the influence of the field- 
induced "effective masses" of the particles on kinematics. 
It is shown that a novel channel of the photon decay 
7 — >■ 7+a is opened kinematically. In Sec. IV the effective 
077 vertex is used to study the axion decay into two 
photons a — > 7 + 7. A catalyzing influence of the external 
field on this process is analyzed. In Sec. V we study the 
photon decay 7 — > 7 + a. The energy loss by a photon gas 
due to this process is investigated as an additonal source 
of the axion luminosity in a supernova explosion. In the 
appendixes, the important details of the calculation of 
the 077 vertex are collected. 

II. a 77 VERTEX IN THE CROSSED FIELD 

Invisible axion models are classified into two types, de- 
pendingon whether or not axions couple to leptons (see, 
e.g., ]5|,|lT[|). Here we investigate the field-induced 077 in- 
teraction of "Dine-Fischler-Srednicki-Zhitnitsky (DFSZ) 
axions" |^3| which couple with both quarks and leptons 
at the tree level. 

The Lagrangian describing the axion-fermion interac- 
tion is 

£a/ = -ig a f(f75f)a, (3) 




FIG. 1. The effective vertex of the interaction of a pseu- 
doscalar particle with two photons. Double lines correspond 
to fermion propagators in an external electromagnetic field. 

where g a f = CfUif/ f a is a dimensionless Yukawa cou- 
pling constant, Cf is a model-dependent factor |Tj] ], nif 
is the fermion's mass, 75 is the Dirac 7 matrix, and a and 
/ are the axion and fermion fields. 

The two-photon axion interaction is investigated in a 
general case when the photons are not assumed to be 
on the mass shell. It means that the expression for the 
077 vertex can be used as an effective Lagrangian of 
axion-photon interaction in studies of processes involv- 
ing axion or other pseudoscalar particles with a coupling 
of type (||) . In the third order of the perturbation the- 
ory the matrix element of the effective 077 interaction 
in an external field is described by two three-point loop 
diagrams (Fig. [I]), where double lines imply that the in- 
fluence of the external field in the propagators of virtual 
charged fermions is taken into account exactly. 

The general form for the matrix element corresponding 
to the diagrams in Fig. [l] is 

_ 4iraQ 2 f g af 

~ ^2E a V x 2ujiV x 2uj 2 V 

x J d 4 xi d 4 A d^Y exp [i(xx(q\ + q% — p))] 

xe^(-i + fey) + ^(XFY)}) (4) 

x S P {S(Y)(e 2l )S(X - Y)(s l7 )S(-X) 75 } 

+ <-> £2,92), 

where a = 1/137 is the fine-structure constant; e > is 
the elementary charge, Q f is a relative fermion charge in 
the loop; p is the four-momentum of the decaying axion, 
9i,2 = (w, q)i,2 and ei t % are the four-momenta of the final 
photons and their polarization four-vectors, respectively; 
(£7) = £^7^, 7 P are the Dirac 7 matrices; X — x\ — x 2 , 
Y = X1 — X3; S(X) is the translationally invariant part of 
the fermion propagator in the crossed field (Appendix |A|) . 

Further we will calculate a vertex function A connected 
with the matrix element Eq. (Q) in the following way: 

„ WJWQ, -„-.,) , 

^2E a V x 2wi^ x 2uj 2 V K ' 

The function A, thus defined, is the effective Lagrangian 
of the axion-photon interaction in the momentum repre- 
sentation. The calculation technique of the three-point 
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vertex 077 in the crossed field is given in Appendixes [B| 
and ^|. In the following, we will consider only the field- 
induced contribution A^ F ) = A — A^°\ where A^ ) is the 
vacuum part of A, to the effective Lagrangian, corre- 
sponding to the transition a(p) — > 7(91) + 7(92): 



A (F) = « 

f 



Q}9af 



(/l/ 2 )J + (/l^/ 2 )Jl 



+ (/l^)(/2^)J2 + (/2^)(/l^)J 2 



(6) 



(qiTTq2) {qiTTqi) 



T - ElilL T - - r T F, - m f 

•> i-iv p i •> [tv — 2 tlvctpJ a/3 i -r/ — g Q j 



(7) 



Jo = i J dz J dxdy [cos ip2 

D 

00 

Ji = 2 J zdz J dxdy 

D 

X 



f 



f 



X2 



sin 1^2 e 



J2 = — / z 2 dz / dxdy y(l — 2x — y) cos ^2 e 
X2 J J 
D 

OC 

J zdz J dxdy (I — 2x) sin^e"'*, 
i = 2, 3, 



J 3 = 



(8) 



Ji = J%{q\ —> 92), 
<P = ztpo + —ip l7 



ip = l- 2xy 



(9192) 
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2/(1 -y)—2-x(i-%)—o 



2 • 



fi = - x)\xi\ - 2/(1 - y)|xi|] 
+ 2xy(I - x - y) 2 (|xiX 2 | - X1X2) 
+ 2xy[xy + (x + y)(l -x- y)]{\xiX2\ + X1X2) > 0, 



tp 2 = -2z 2 xy(l -x-y) 
2 _ (pFFp) 2 _ 



{qiTTqi) 



1,2. 



nr. 



The integration area D in the integrals Ji is x, y > 0, 
£ + V < 1- 

We note that the vertex A^ F ) (|) contains as a limit 
an amplitude of the field-induced effective 07 interaction 
which can be used in studies of axion-photon oscillations. 
The amplitude of the a(p) —> 7(g) transition (the dia- 
gram is shown in Fig. is obtained from Eq. (||) by 



a(p) 




FIG. 2. The diagram describing the effective axion-photon 
interaction. 



means of replacing the field tensor of one of the photons 
by the electromagnetic field tensor: 

qia — > q a , fla/3 — > /a/3, 92a — > 0, f^ a j3 — » —lF a p. 

As a result, the amplitude of the a — > 7 transition can 
be presented in the form []: 



M = -^(/F)_ 

Z7T 771/ 



0)], (9) 



7 = 



2/3 



ir/2 



/fa) sin" 1 / 3 ^, 



/(?/) = i du exp < — i 7771 + 



v= 1- 



9 ■ 2 
<7 sin 1 

Am 2 



XSin 



2/3 



where /(77J is the Hardy-Stokes function and v is the 
dynamic parameter of the axion defined in Eq. (j8|) . 



III. FIELD-INDUCED "EFFECTIVE MASSES" 
OF THE PARTICLES 

The vertex 077 Eq. (||) can be used as an effective 
axion-photon Lagrangian in studies of processes such as 
the radiative decay a — > 77, photon splitting 7 — > 7a, 
and the Primakoff-type process for photoproduction of 
pseudoscalars on electrons 7 + e — > e + a. 

We note that in a homogeneous electromagnetic field 
the energy-momentum conservation law for axion-photon 
processes coincides formally with the vacuum one 
However, in calculating the probabilities of the processes 
one has to integrate over the phase space of the final par- 
ticles taking into account their nontrivial field-induced 
kinematics. This is due to the fact that an external field 



^For a numerical evaluation note that the external field ten- 
sor F^v is presented here in rationalized natural units. 
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FIG. 3. The photon dispersion curves in the crossed field. 



FIG. 4. The diagram of the axion transition a 
external electromagnetic field. 



has a negative minimum and changes its sign with in- 
creasing the dynamic parameter \- 

The fact that both of the photon eigenmodes have neg- 
ative "squared masses" in the region of x 15—18 means 
that two-photon decay is kinematicahy opened even for 
massless pseudoscalars (e.g., arions The difference 
in values of the field-induced "squared masses" of the 
first and second photon eigenmodes makes possible the 
process of the photon splitting 7 — > 7a, where a is an 
arbitrary relatively light pseudoscalar. 

In general, it is necessary to analyze the influence of the 
external field on the axion mass m a . The field-induced 
contribution to m a is connected with the real part of the 
amplitude AM of a — > // — > a transition via the fermion 
loop by the relation: 



5ml 



-Re AM. 



(12) 



plays the role of a peculiar medium with dispersion and 
absorption. In our case the external crossed field plays 
the role of a homogeneous anisotropic medium. The pho- 
ton "effective masses" squared fi\ induced by the external 
field are defined as the eigenvalues of the photon polar- 
ization operator: 



n„ 



3 

E 

A=l 



b w b w 



W^fAA'fW. (10) 



To obtain a correct result for AM one has to use the 
Lagrangian with the derivative as it was first emphasized 
by Raffelt and Seckel Epl: 



(13) 



In the second order of the perturbation theory the am- 
plitude a — > // — > a is described by the diagram in Fig. ^ 
and can be presented in the form: 



Here are the polarization operator eigenvectors de- 
scribed in Appendix [5]. We stress that only two eigen- 
modes of the photon propagation with polarization vec- 
tors 



b {1) 



~ ^(qFFq) ' 
(qF) a 
(qFFq) 



(11) 



are realized in an external electromagnetic field. The 
analysis of the photon polarization operator H uv in a 



one-loop approximation in the crossed field 24 shows 
that the dispersion curves corresponding to the above- 
mentioned photon eigenmodes (Fig. ||), though being 
alike in their qualitative behavior, are different quanti- 
tatively. It is seen from Fig. |3| that each dispersion curve 



AM 



(14) 



xSp[S(-Z) (qj)KS(Z) (97)75], 



where Z = x — y. The integration with respect to the 
four-coordinate Z is r educ ed to the generalized Gaussian 
integrals of the type ( |C6|) from Appendix |(] and results 
in 



AM 



16tt 2 



7T/2 



2q 2 



7T/2 



sin ( 



dt 

T 



2 2 

m f x 



$ = $0 + 



t 3 x 2 sin 4 1 
48 



t dt 



$n = t 1 - 



(15) 



o . 2 

q sin 1 
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The main contribution to Srr? a has the following asymp- 
totic behavior at large values of the dynamic parameter: 



^-|g|r 2 (2/3)(6 X ) 2 / 3 , 



(16) 



and permits one to estimate the field-induced correction 
to the axion mass: 



<5m 2 



10" 



' C f mj 
lMeV 2 



2/3 



(17) 



Here we used that m a and g a f are connected with the PQ 
symmetry-breaking scale f a by virtue of the relationships 
m„ = 0.62eV x (10 7 GeV// a ) and g af = C f m f /f a @. 
Equation ( |l7j ) shows that we can neglect the influence of 
an external electromagnetic field on m a . 



IV. THE AXION DECAY a -> 7 + 7 

A. The decay amplitude 

The expression for the effective 077 interaction (^J) is 
used to study the ultrarelativistic (E a 3> m a ) axion de- 
cay a(p) — > j(qi) + 7(92)- Because of the smallness of 
the axion mass (m a < 10 -2 eV), even axions with the 
energy of several eV are ultrarelativistic. Assuming that 
the field-induced photon "masses" are also small in com- 
parison with their energies, below we will consider pho- 
tons as massless. This case corresponds to the collinear 
kinematics: 

Before starting to study the physically most interesting 
case of the ultrarelativistic axion decay let us analyze all 
tensor structures of the effective axion-photon vertex (||), 
which define the decay amplitude, in an attempt to reveal 
the dominating contribution. Note that for real photons 
the terms with J 3 and J3 are suppressed by the smallness 
of the field-induced photon "masses" (<?i/i) Q ~ ^ -Sa- 
lt is natural to denote the other terms in Eq. (|) asM(°), 
M' 1 ), and M^ 2 ) depending on the power of the external 
field tensor. 

Due to the Lorentz invariance the amplitude can be an- 
alyzed in any frame. We will analyze it in the rest frame 
of the decaying axion with the components of the four- 
momenta p, q±, and q 2 being of order m a . In this case 
it is sufficient to allow for the order of the dimensional 
quantities M (0 \ and M (2 ^. 

m f m f 



2 
— 



(18) 



~ -\(f 1 F)(f 2 F) ~ F 1 
f f 



where F 1 stands for the strengths of the magnetic and 
electric fields in the axion rest frame. In this case the elec- 
tromagnetic field in Eq. ([IJ|) is obtained by the Lorentz 
transformation from the laboratory frame, in which the 
external field F is given, to the rest frame of the decaying 



m a 



In view of Eq. ( |l9j ) , the expressions 
in the form 



(19) 



181) can be written 



m a E a 



in 



F, 



(20) 



/ 



p2 



The following is seen from Eqs. (|20|). 

1. The external electromagnetic field affects substan- 
tially the ultrarelativistic axion decay because the 
field reduces the suppression caused by the small- 
ness of the axion mass. A similar catalyzing effect 
of the external electromagnetic fields of various con- 
figurations on the radiative decay of the massive 
neutrino was discovered in |2q,|27 



2. The contribution to the amplitude bilinear in the 
external field dominates because it does not contain 
the above-mentioned suppression factor. 

With regard to the analysis performed only terms bilinear 
in the external field should be kept in the amplitude of 
the ultrarelativistic axion decay: 



M 



E 



Qj9af 



(21) 



\ < —(fiF)(f 2 F)J+(e liqi «-> £ 2 ,g 2 ) 

.X2 



1 l-x 

J(Xi,X2) = Jdx J dyy(l-y-2x)i] 3 (l-7]f(7])), 


V = { W ~ x)X2 - 2/(1 - v)xif (22) 

1 ~ 1/3 

+ ixy [xy + (1 - x - y)(x + y)] X1X2 J 

where f(f}) is the Hardy-Stokes function defined in 
Eq- @t Xj Xii an< i X2 are the dynamic parameters de- 
termined in Eq. (||). We have neglected terms of order 
m 2 /£" 2 , fif/E 2 in the argument of the function f(rf). 
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B. The decay probability 

The probability of the field-induced decay of the axion 
has the following form: 



1 1 



d 3 gid 3 q 2 
2! 32ir 2 E a J wiw 2 



S^(p- qi -q 2 )J2\M\ 



A, A' 



32irE, 



— [ dtY \M 



(23) 



The factor 1/2! in the integration over the final state 
takes into account the identity of two final photons. The 
process considered is two-body decay, but its amplitude 
is not a constant as in vacuum because it depends on 
the external electromagnetic field tensor. That is why 
the calculation of the decay probability is not reduced 
to multiplying the matrix element squared with the final 
phase space. Under the summation of the polarizations 
A, A' in Eq. (^3|) it is convenient to use the polarization 
vectors The expression for the probability is 



1 /a\ 2 



dt r 



(24) 



Y,Q 2 f9a f m fX 2 J(t X ,(l-t)x) 
f 



With the strong hierarchy of the fermion masses, 

6 



A 
4 



the contribution from the fermion with the maximum 
value of the dynamic parameter \ 2 = e / (pEFp) /m^ can 



dominate in Eq. ( |24[ ) 

WW 



3.32 



a\ 2 (Qjdafmf) 2 



nE a 



P(X)- 



(25) 



The plot of P(x) is shown in Fig. 0. We also give the 
asymptotic behavior of the function P(x) for both small 
and large values of the dynamic parameter x : 



P(X) 

P(x) 



2.31 x 10 _5 x 8 (l + 9.5x 



.8 32.0 



x i/ 3>1 



X ^ X^ 



As it is seen from the asymptotic behavior of the func- 
tion P(x) at large x, the decay probability (^Bj) does not 
depend on the dynamic parameter. Note that a similar 
behavior of the amplitude and probability at x 3> 1 takes 
place in the photon splitting 7 — > 77 p9| . 




FIG. 5. The plot of the decay factor P(x) as a function of 
the dynamic parameter \- 



To illustrate the catalyzing influence of the external 
field on the ultrarelativistic axion decay, let us compare 
Eq. (|25| ) with the well-known axion decay probability in 
vacuum |7 11 1: 



W n = 



2 4 

6^E a ' 



(26) 



Here g ai = (a/27r/ a ) x (E/N - 1.92 ± 0.08) Q, where 
E and TV are model-dependent coefficients of the electro- 
magnetic and color anomalies. A comparison of Eqs. (I2T 
and ||), 



R 



f 



Wn 



2.12 x 10^ 



P(X), (27) 



demonstrates a strong catalyzing influence of the exter- 
nal field on the ultrarelativistic axion decay (E a ^> m a ), 
because has no suppression factor associated with 

the smallness of the axion mass, except the coupling con- 
stants (g a f,ga-y ~ fa 1 )- Let us estimate in Eq. ( p7j ) the 
contribution from the electron in the loop because this 
fermion is the most sensitive to the external field: 



Re ~ 10 



37 / COS 2 P 



\E/N - 1.92 



10- J eV 



^(X)- (28) 



Here cos j3 determines the electron Yukawa coupling 
g ae = 3 cos 2 /3 (m e /f a ) in the DFSZ model |^3| . Large 
values of the dynamic parameter X 1, when P(x) ~ lj 
can be realized, for example, in the case of the decaying 
axion energy ~ 10 MeV and the magnetic-field strength 
B - B e ~ 4.41 x 10 13 G. 



V. PHOTON SPLITTING 7 -» 7 + a 

As it was pointed in Sec. Ill the external field has a 
substantial influence on the photon dispersion relation so 
that the photon splitting 7 — > 7a is allowed kinematically 
due to the difference in the field-induced effective photon 
"masses" . This process can play a role as a possible addi- 
tional mechanism of energy loss by astrophysical objects. 



G 



As this process is a crossing channel of the axion decay 
a (p) ~ * 7(91) + 7(92) and its amplitude can be easily 
obtained from the amplitude ( plf ) using the replacement: 



Qi -* -<?i> p->-p, 
which corresponds to 



/1 



-/1 



-Xi- 



(29) 



(30) 



The kinematics of this process is close to the collinear one 
as well as in the ultrarelativistic axion decay. Note that 
in studying this process in different ranges of the dynamic 
parameter one should take into account the specific dis- 
persion behavior of the photon modes in the initial and 
final states. 



A. The limit of small values of the dynamic 
parameter xi *C 1 

In the case of small values of the initial photon dynamic 
parameter xi the splitting of the photon with the first 
polarization J] is allowed kinematically only due to the 
condition [if > /i?, (see Fig. ||): 
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(i) 



(31) 



With the photon polarization vectors ( |ll| ) the amplitude 
can be presented in the form: 



4a . 



: —t (l-t)J2 Q)9a f m f x\ J(tXl,Xi), (32) 
/ 



where t = a^/^i is the relative energy of the final pho- 
ton. The function J defined in Eq. ( p2| ) has the following 
asymptotic behavior at small values of its arguments: 



^(*Xi>Xi) 



-Xi (l-2t) + 0(xi). 



Xi<l 

The probability of the photon splitting is 
1 

dt\M\ 2 



1 



167TCJ1 



(33) 



4.8 x 10-° ( - 

. 7T 



TTLOl 



To illustrate a possible application of the result ob- 
tained we estimate the contribution of this process to 
the axion emissivity Q a of the photon gas: 



* Numbering of the photon types is made by their polariza- 
tions in accordance with Eq. (p"i"|). 



d 3 q 



(2tt) 3 



(34) 



dt—— (l-t) [1 + n B {uit)], 
dt 



where ns(^>i) and ns(uJit) are the Planck distribution 
functions of the initial and final photons at temperature 
T, respectively. In Eq. ([m]) we have taken into account 
that the photon of only one polarization (the first one in 
this case) splits. Note that the dynamic parameter xi 
in dW( p > I dt depends on the photon energy u)\ and the 
angle 9 between the initial photon's momentum qi and 
the magnetic- field strength B: 

wi B 
xi = — -5- sm e > 

mj B f 

where Bf = m 2 / eQ / is the critical value of the magnetic- 
field strength for the given fermion. Neglecting the pho- 
ton "effective masses" squared (uj 2 — |qi| 2 + £i 2 — |qi| 2 ), 
the result of the calculation of the axion emissivity ( |34| ) 
can be presented in the form 




T n B 8 . 



(35) 



The numerical factor C is expressed in terms of an infinite 
sum of the Hurwitz £ functions £(fc, z) = J] + z ) ■ 

n=0 



C = 



63 2 



sin y 9d8 / dtt z (1-f) (l-2t) z 



x 10 dx 



(e x - 1)(1 - e' xt ) 


2 13 ~/C(3,n) 4£(5,n) 3C(7,n) \ 
= w L {—s TrfT- + J = 2 ' 15145 ' 

n— 1 v 7 

where x = u>i/T is the relative energy of the decaying 
photon. 

Let us estimate the contribution of the photon splitting 
7 — > 7a into the axion luminosity in a supernova explo- 
sion from a region of order of a hundred kilometers in 
size outside the neutrinosphere. In this region a rather 
rarefied plasma with the temperature of order of MeV 
and the magnetic field of order 10 13 G can exist. Under 
these conditions the axion luminosity is 



^ 3xlo 7lio- 



T 



lMeV 



(36) 



B 



10 13 G 



R 



10 3 km 



7 



The comparison of Eq. (|36]) with the total neutrino lu- 
minosity L v ~ 10 52 erg/s from the neutrinosphere shows 
that the contribution of the photon splitting to the en- 
ergy loss is negligibly small. Here we used the strongest 
restriction on the axion coupling constant with electrons 
is ffae ~ 10- 13 IPH- 



B. The limit of large values of the dynamic 
parameter X\~^> 1 

At large values of the dynamic parameter xi the field- 
induced "effective mass" squared of the photon of the sec- 
ond polarization becomes greater than that of the photon 
of the first polarization (/x 2 . > fif, see Fig. ||). In this case 
the channel 



7 



(2) 



7 



(i) 



(37) 



is kinematically opened only. The process amplitude at 
large \i is 



M~ — (l-t)J2 Qf9afm f I(t), 
n f 

y (1 — y — 2x) dxdy 



(38) 



I(t) = 



D 



[x(l - x)t - 2/(1 - y)] 2 + Axyt (1 - x - y) 2 



In the third order of the perturbation theory, the field- 
induced vertex 077 (^) was obtained for the case when 
all the particles were off the mass shell. This vertex can 
be used as the effective Lagrangian of the 077 interaction 
in the investigation of processes involving axion or other 
pseudoscalar particles with a coupling of the type (||) . 

As a specific example of the way this vertex can be 
used we considered the axion decay a — > 7 + 7 in an 
external field. The dominant contribution to this pro- 
cess (not suppressed by the axion mass) is presented in 
Eq. ( pl| ) . The decay probability ( p5|) is calculated in the 
limits of both large and small values of the dynamical pa- 
rameter x- Comparison of Eq. (^) with the probability 
of the two-photon axion decay in vacuum shows a strong 
effect of the external field (~ 10 37 ) on the decay of the 
ultrarelativistic axion (E a ^> m a ). 

As another example, we have studied the photon split- 
ting 7 — > 7 + a which could be of interest as an addi- 
tional mechanism of energy losses by astrophysical ob- 
jects. This forbidden in vacuum process becomes kine- 
matically possible because photons of different polar- 
izations obtain different field-induced effective "squared 
masses" . At the same time an external field influence on 
the axion mass is negligible. The axion luminosity esti- 
mated for supernova conditions turn out to be small in 
comparison with the neutrino luminosity. 



where the integration is carried out over the area D: 
x,y > 0, x + y < 1. It is seen from Eq. ( |38"| ) that the 
amplitude does not depend on the dynamic parameter. 

With the amplitude ([38j) the splitting probability has 
the form: 



f m f 



ttuji 



(39) 



In real astrophysical objects with strong magnetic 
fields the temperature T > 10 MeV can exist in the cen- 
tral regions where from both components of the active 
medium, the magnetic field and a plasma, the plasma 
component substantially dominates. Thus, the expres- 
sions obtained are not applied for the estimation of the 
luminosity of astrophysical objects because we have taken 
into account the influence of the external field only. 



VI. CONCLUSION 

In this work we have investigated the effective 077 in- 
teraction (a is a pseudoscalar particle) induced by an 
external constant crossed field. Calculations of processes 
in this electromagnetic field configuration are very gen- 
eral, being a relativistic limit of calculations in an arbi- 
trary weak smooth field. For the pseudoscalar particle 
wc considered the axion, the most widely discussed par- 
ticle corresponding to the spontaneous breaking of the 
PQ symmetry. 
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APPENDIX A: THE PROPAGATOR OF A 
CHARGED FERMION IN A CROSSED FIELD 

A charged fermion propagator S^ F \x,y) in the exter- 
nal crossed field in the proper time formalism [ p9| , j30[ c an 
be presented in the relativistically invariant form p7| 



S (F) (x,y) 
S(X) 



= J^{x,y) 



-e f I d£ 



S(X), 



(Al) 



M0 + 7;F^~y) 1 



167T^ J ,s 





5- (^7) - ^(^7)75 

Is 2 



se t 



-J-(XFF7)+m f 



sm f e f 



( 7 j? 7 ) 



8 




FIG. 6. n-point loop fermion diagram with cut legs in an 
external electromagnetic field. 



x exp 



sm. 



X*_ 

4s 



sej 
12 



(XFFX) 



where — (x — y)^; rrif and e/ are the mass and charge 
of the virtual fermion (e t = eQ f , e > is the elementary 
charge, Qf is the relative fermion charge); A^, F^, and 
F^ u are the four-potential, tensor, and dual tensor of 
the crossed field; (Xj) — A M 7 M , 7 M , 75 are the Dirac 7 
matrices. 

The inclusion of translationally noninvariant factors 
f2(xi, cCi+i) arising from the corresponding propagators 
S^ F \xi, Xi+i) in calculations of n-point loop diagrams is 
discussed in detail in Appendix pi 



APPENDIX B: INCLUSION OF 
TRANSLATIONALLY NONINVARIANT 
FACTORS Q(x,y) OF FERMION PROPAGATORS 
IN n-POINT LOOP DIAGRAM 



In this Appendix we discuss the inclusion of the trans- 
lationally noninvar iant factor Q(x,y) from the propaga- 
tor S^ F \x, y) Eq. (Al) in the calculation of n-point loop 
diagrams. In particular, the diagrams in Fig. |l] describ- 
ing the effective 077 interaction are three-point loop di- 
agrams. 

The four-potential A^(x) of a constant external elec- 
tromagnetic field can be written as 



.4 



l ( - 1 " 

I' ~ 7^\ x F)i_i — ~^X Fufj,. 



(Bl) 



With Eq. (Bl), the expression for £l(xi,x 2 ) is substan- 
tially simplified: 



Q{x x ,x 2 ) = -e } / 



M0 + \f^-x2)" 



e -^(xiFx 2 ) 



(B2) 



The contribution from t he translationally noninvariant 
factors £l(xi, Xi+x) Eq. (B2) to the amplitude of an n- 
point diagram in Fig. ^ is generally given by 



fit,* = ~k y^xxtFx l+1 ) 



n-l 



J2(Z l FZ l+1 ) 



.r„ + i=xi 



i = l 



X% 2-2 + 1, 



Z n = x n -xi, Zj = 0. (B3) 



We note that, while for each propagator S^(x,i,Xi + i) 
the phase fl(xi, ajj+i) is generally translationally and 
gauge noninvariant, the total phase of the n-point fer- 
mion loop Qtot Eq. (B3) only depends on the difference 
of the coordinates and satisfies the requirement of gauge 
and translational invariance. 

The simplest situation is realized for the two-point loop 
diagram when the total phase is zero: 



Qtot = fl(xi,x 2 ) + Q(x2,Xi) 
X = x\ — x 2 . 



(XFX) = 0, 



For the three-point loop diagram studied in this paper 
we have 

fit* = fl{ Xl ,x 3 ) + n(x 3 ,x 2 ) + n(x 2 ,xi) = -^-(XFY), 
X = x\ — x 2 , Y — x\ — X3. 



APPENDIX C: CALCULATION OF THE 
VERTEX FUNCTION A 

In this Appendix we give the basic stages of the calcu- 
lation of the three-point vertex function A describing the 
effective 077 interaction (a is an arbitrary pseudoscalar 
particle) in the external crossed field. After integration 
with respect to the coordinate X\ in Eq. ^ and separa- 
tion of the delta- function S^(p — qi —q 2 ), the expression 
for A can be represented according to Eq. (0) as 

A = -ie)g af J d 4 X d 4 Y (CI) 

x exp (-» [( gi X) + (q 2 Y) + ^1{XFY))) 

x Sp{S(Y)(e 27 )S(X - y)( £l7 )S(-A)7 5 } 
+ {sx,qi <-> s 2 ,q 2 ). 

With the explicit expression for the propagators ( |Al[ ) the 
vertex function is 

2 00 

. _ a QjOafny f dsdvdr -i m 2 f (s + v + r) 
~ it ^ 8(4tt) 4 J s 2 v 2 t 2 6 



x J d 4 Xd 4 Y I ((A - Y)R 1 X) + (YR 2 (Y - A)) 

™2 



111 f 

(YR 3 X) - -J-R 



(C2) 



X 2 (\ 1 



x exp ( - i\ {qiX) + (q 2 Y) + _!_ + _ 
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y2 

4 

„2 



1 1 



(XY) 



2 y 



^(s + v) (YFFY) 
(ei,gi <-> £2,92), 



(XFFY) - — (XFY) 



where s, v and r are proper time variables defined by 
three propagators of the form (Al), X — x\ — x 2 , 
In Eq. (F~ 



Y = x\ — X3. In Eq. (|C2j) all the dependence on four- 
coordinates X and Y is shown explicitly in the preexpo- 
nential factor and in the exponent, so that Rq and the 
tensors (Ri) a p {i = 1,2,3) are functions of proper time 
variables and the constant external field tensor: 



Rq 

(Rx)afd 
{R 2 )aj3 
{R?,)al3 

hx(s) 
h 2a {s) 



Sp{75Ms) ( £ 27) hi{v) (£17) hi(r)}, 
Sp{75M s ) ( £ 27) h 2a (v) (£17) h2p(r)}, 
Sp{ 75/11 (t) h 2a (s) (£27) h 2 [j{v) (£17)}, 
Sp{ 75/11 («) (£17) ^(t) /i 2q (s) (£27)}: 
2 + e /S ( 7 F 7 ), 

1 



2,s 



7a 



(lFF) a 



The integration with respect to the four-coordinates X 
and Y in Eq. (C2) reduces to the calculation of the gen- 
eralized Gaussian integrals of three types: scalar, vector, 
and tensor of rank two. The calculation method of such 
integrals which arise in the inte gra tion of the two-point 
loop diagrams is given in Ref. p7[ . However the inte- 
gration with respect to the coordinates is significantly 
simplified by formally introducing the eight-coordinate Z 
and the eight-momentum Q defined in terms of the four- 
coordinates X, Y and the four-momenta of photons q%, 
q 2 in the following way: 



Y„ 



Ql 



9V 



(C3) 



In terms of the generalized coordinate Z^ and the gen- 
eralized momentum the vertex function ( p2| ) can be 
written in the most compact form: 



2 

, _ a y> Qf9afmf f dsdvdr _i TO 2( s + w + r ) 

~ 7T 2f 8(47r) 4 J S 2 V 2 T 2 6 

R + (ZNZ) 



x J d s Ze -i[\(ZGZ) + (QZ)] 
+ (ei,qi «-> s 2 ,q 2 ). 



(C4) 



The matrices G and N can be presented as the Kronecker 
product of (2 x 2) matrices and Lorentz tensors of rank 
two [(4 x 4)-matrices] : 

G aP = g aP Gx +e f F a[3 G 2 + e) (FF) a p G 3 , (C5) 

q _ 1 / s(v + T) —ST 

1 ~~ svt V -ST t(s + v) 



G 2 = 



1 

-1 



, -, 1 / V + T -V 

-„ S + V 



N aP = Nx(Rx) a f) + N 2 {R 2 ) aP + N 3 (R 3 ) af3 , 



1 

-1 



No 





-1 1 



, N 3 





1 



Thus, t he i ntegration with respect to the coordinate Z^ 
in Eq. (C4) has been reduced to the calculation of the 
generalized Gaussian integrals of two types - the scalar 
and tensor of rank two: 



J= [ d s Ze -i [\(ZGZ) + {QZ)] 



(C6) 



-(47T) 4 



AQG-'Q) 



o 2 j 

dQ^dQ" 



(s + v + t) 2 
T,^Jd^ZZ^e-ni(ZGZ) + (QZ)] = 

= [4(QG- 1 ) tl (G- 1 Q) v -2iG-Z] J, 
where G _1 is the inverse of the matrix (|C5|) : 

G~p = g aP Gx +e f F aP G 2 + e) {FF) a p G 3 , (C7) 



_ 1 f T(S + V) ST \ f, _ svt ( -1 

Gl -~z{ st s (v + t) h G2 - — {i r 



(G3)n = — ~ z(4sv — s 2 — v 2 ) — 6svt 



3z 2 

(^3)12 = (G3J21 



3z 2 



V 3 + S 3 + T 3 + 3STV - Z(s 2 + T 2 ) 



(G 3 ) 2 2 = z(Avt - t 2 - v 2 ) - 6SVT 



3z< 



Here the variable z is defined as z = 



s + v + T. 



APPENDIX D: EXPANSION IN THE BASIS IN 
AN EXTERNAL ELECTROMAGNETIC FIELD 

Calculations are significantly simplified by using a ba- 
sis which is constructed of the photon four-momentum 
and the external electromagnetic field tensor F^ Q : 



&« = (gF) M , bf = (qF)^ (Dl) 
bf = q 2 {qFF)^ - q^qFFq), = q„. 



In the specified basis a polarization vector of an arbitrary 
photon is 



(Ff) 
2(qFFq) 



(Fq\ 



(Ff) 
2(qFFq) 



(Fq)» (D2) 



(qfFFq) 
q 2 (qFFq) 2 



q 2 {FFq) ll -{qFFq)q ll 
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where f a /3 = q a £p — qp^a i s the photon field tensor. Note 
that the last term in Eq. (D2) is not equal to zero only 
for virtual photons. 

Taking into account Eq. (|D2] ) , we can reduce all struc- 
tures of the type (hfl 2 ) that appear in the expression for 
the vertex function A to the form: 



5] R.D. Peccei, Report No. hep-ph/9606475 (unpublished) 



(hfk) = (fF) 



(hq)(qFl 2 ) - (l 2 q)(qFh) 
2{qFFq) 

(hq){qFh) 



(D3) 



2(qFFq) 
(hqKqfFFh) - (hg)(qfFFl 2 ) 
(qFFq) 

where l± and l 2 are arbitrary four-vectors. Below we 
give som e structures for which we used the representa- 
tion (IfF 



( gi f 2 FF qi ) 



(qiFFq 2 ) 
2(q 2 FFq 2 ) 

(f 2 F)( qi Fq 2 ) + (f 2 F)( qi Fq 2 ) 



(D4) 



(qif 2 FFq 2 ) = - [(f 2 F)( qi Fq 2 ) + (f 2 F)( qi Fq 2 ) 
(qif2Fq 2 ) = - (qiq 2 )(f 2 F), 



where (fi) a /3 is the electromagnetic field tensor of the ith 
photon (i — 1,2). 

Note that the following relationship: 



mf 2 )^ + (T 2 f x ) 



(D5) 



proves very helpful between two arbitrary antisymmet- 
ric tensors (Ti) M „ and [T 2 )^ v . If for tensors T\ and 
T 2 we take the constant electromagnetic field tensor 
T\ = T 2 = F, which meets the condition (FF) = 0, 
then from Eq. (D5) the condition on the electromagnetic 
field tensor F follows automatically, 

(FF) MV = 0. 
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